Hierarchical Mean-Field Theories in Quantum Statistical Mechanics 
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We present a theoretical framework and a calculational scheme to study the coexistence and competition of 
thermodynamic phases in quantum statistical mechanics. The crux of the method is the realization that the 
microscopic Hamiltonian, modeling the system, can always be written in a hierarchical operator language that 
unveils all symmetry generators of the problem and, thus, possible thermodynamic phases. In general one 
cannot compute the thermodynamic or zero-temperature properties exactly and an approximate scheme named 
"hierarchical mean-field approach" is introduced. This approach treats all possible competing orders on an equal 
footing. We illustrate the methodology by determining the phase diagram and quantum critical point of a bosonic 
lattice model which displays coexistence and competition between antiferromagnetism and superfluidity. 
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PACS numbers: 05.30.-d, 05.70.Fh, 05.30.Jp, 75.10.Jm 

During last decade we have witnessed great advances in 
materials synthesis and fabrication. The rate at which new 
compounds with multiplicity of distinct phases and character- 
istic functionalities are generated has outpaced our complete 
physical understanding of the fundamental principles behind 
such complex phenomena. For example, whether the mecha- 
nism controlling the coexistence and/or competition between 
magnetism and superconductivity or Bose-Einstein condensa- 
tion has the same physical origin in different classes of ma- 
terials is still an open question [|l]]. On the other hand, the 
possibility of control and tunability of the interactions of the 
elementary constituents (i.e, quantum control) offers the po- 
tential to design new states of matter with unforeseen appli- 
cations [^]. Despite great theoretical advances there is a lack 
of a systematic and reliable methodology to study and predict 
the behavior of these complex systems. It is the main purpose 
of this paper to present a promising step in that direction. 

The modern theory of phase transitions starts with Landau's 
pioneering work in 1937 [|^. One of his achievements was 
the realization of the fundamental relation between sponta- 
neous symmetry breaking and the order parameter (OP) that 
measures this violation, thus giving simple prescriptions to 
describe order in terms of irreducible representations of the 
symmetry group involved. Another was the development of a 
phenomenological calculational scheme to study the behavior 
of systems near a phase transition. Landau's theory has been 
successfully applied to study phase transitions where thermal 
fluctuations are most relevant. Certainly, the theory was not 
designed to study zero-temperature (quantum) phase transi- 
tions, i.e., the qualitative changes of the macroscopic state of 
the system induced by tuning parameters of its Hamiltonian. 

In the quantum description of matter, a physical system is 
naturally associated with a language of operators [^. In pre- 
vious work [|| ^ we outlined a framework to identify OPs 
based upon isomorphic mappings to a hierarchical language 
(HL) defined by the set of operators which in the fundamen- 
tal representation (of dimension D) has the largest number of 
symmetry generators of the group. Any local operator can be 
expressed as a linear combination of the generators of the HL. 
The building of the HL depends upon the dimension D of the 
local Hilbert space modeling the physical phenomena one is 
investigating. For instance, if one is modeling a doped anti- 



ferromagnetic (AF) insulator with a t-J Hamiltonian [^, then 
D=3 (i.e., there are three possible states per site) and a HL is 
generated by a basis of su{3) in the fundamental representa- 
tion [||, ^]. As explained and proved in Refs. |]], there is 
always a HL associated to each physical problem. These ideas 
complement Landau's concept of an OP providing a mecha- 
nism to reveal them, something that is outside the groundwork 
of Landau's theory. Indeed, this theory does not say what the 
OPs should be in a general situation. 

It turns out that these isomorphic mappings not only un- 
veil hidden symmetries of the original physical system but 
also manifestly establish equivalences between seemingly un- 
related physical phenomena. Nonetheless, this is not sufficient 
to determine the exact phase diagram of the problem: One 
has to resort to either numerical simulations with their well- 
known limitations or, as will be shown in the present paper, 
to a guided approximation which at least preserves the qual- 
itative nature of the possible thermodynamic states. A key 
observation in this regard is the fact that typical model Hamil- 
tonian operators written in the HL become quadratic in the 
symmetry generators of the hierarchical group, and this result 
is independent of the group of symmetries of the Hamiltonian. 

This latter result suggests a simple approximation, based 
upon group theoretical grounds, which deals with competing 
orders on an equal footing and will be termed hierarchical 
mean-field theory (HMFT). In a sense, that will become clear 
below, HMFT constitutes the optimum mean-field (MF) or 
saddle-point solution that approximates the energy and cor- 
relation functions of the original problem. The HMFT is dis- 
tinctly suitable when the various phases displayed by a system 
are the result of competing interactions and non-linear cou- 
plings of their constituents matter fields. From the theoretical 
standpoint these systems are strongly correlated since no ob- 
vious small coupling constant exists, as a consequence they 
exhibit high sensitivity to small parameter changes. It is then 
clear the importance of developing a methodology that treats 
all possible competing orders on an equal footing. 

We will now illustrate the methodology by example and de- 
termine the zero temperature phase diagram of a simple model 
which displays coexistence and competition between antifer- 
romagnetism and Bose-Einstein condensation (superfluidity). 
The model represents a gas of interacting two-flavor (<t =t, |) 
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hard-core bosons with Hamiltonian (t > 0) 



(1) 



where (i, j) stands for nearest-neighbor sites (bond) in an oth- 
erwise regular A^^-sites lattice of coordination z and D = 3. 
The number operator nj ~ Uj-^ + nji (njo- = fej^fej^.), and 



Sj = ^^ln'^^f'^ju is a s = i operator (cr denoting Pauli ma- 
trices). The algebra satisfied by the hard-core bosons is [||]: 
K,biJ = 0, [k.,bU = - 2ni, - 7li^) (if a = a'), 
or — <5ij6j^,6ij^ (if a ^ cr'). Notice that H is an extended 
J-like model of hard-core bosons instead of constrained 
fermions [[Tc|]. These hard-core bosons could represent three- 
state atoms, like the ones used in trapped Bose-Einstein con- 
densates, moving in an optical lattice potential. For the sake 
of clarity we will only consider the AF case J > 0. 

As explained in the introduction the first step in determining 
its phase diagram consists of re-writing iJ in a HL. The latter 
is realized by SU (3)-spin generators in the fundamental rep- 
resentation, and its mapping to the hard-core boson language 
can be compactly written as [0] 
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(2) 



The three components sf = (nji — fijj.)/2, sf = ^ji^jx ^'^^ 



JT ji 

^li^jT generate the spin su{2) subalgebra, i.e., they are 
the components of the local magnetization. The five addi- 
tional components correspond to the Bose-Einstein conden- 
sate and the charge density wave local OPs. In the HL, H 
represents a Heisenberg-like Hamiltonian in the presence 
of an external magnetic field fi' ( J^j^ ~ J^^) 



(3) 
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with Joo ^ V ~ J/2, Joi = Jo2 = t, Jii = Ji2 = J22 = 
J/2, and fi' = §(2^ - J/2) - p.. Note that through the 
mapping we transformed an interacting problem into another 
problem that is quadratic in the basis of the algebra su{3) in 
the fundamental representation, but which is not necessarily 
SU{3) symmetric. This HL furnishes the natural framework 
to analyze the symmetries of the Hamiltonian H. There is al- 
ways an SU{2) spin symmetry generated hy S^^ — S^^, 5^^, 
and S^^. When fi' — and V — 2t, there are five additional 
generators of symmetries related to the charge degrees of free- 
dom. Moreover, if J = V = 2t there is full SU (3) symmetry. 
For /i' ^ 0, the only charge symmetry that remains is a U{1) 
symmetry generated by 5"° (conservation of the total charge). 
In this way the HL, leading to a unique OP from which all 
possible embedded orderings are derived, provides a unified 
description of the possible thermodynamic states of the sys- 
tem. Yet, it remains to establish the orderings that survive as a 



result of tuning the parameters of the Hamiltonian or external 
variables such as temperature and filling. 

For arbitrary values of the parameters J/t, V/t, we do not 
know a priori how to determine exactly the phase diagram of 
H [jn]]. The idea behind the HMFT is to perform an approx- 
imation which deals with all possible local OPs on an equal 
footing with no privileged symmetry axes and, hopefully, re- 
tains the qualitative topology of the phase diagram. With H 
written in the HL one immediately realizes that the simplest 
HMFT can be achieved if we re-write H in terms of SU (3) 
Schwinger-Wigner (SW) bosons (3 flavors a =i,0, t) PI- 
The mapping is expressed as 



/ 71jO - 



Sij) 



with the SW bosons foj^ satisfying the constraint Uji + njo + 
nji = 1. The resulting Hamiltonian (V = 2t with no loss of 
generality) is |g| 



(4) 
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(5) 



where fj, ^ zt — fi and the ordering operators 




which transform as singlets with respect to the genera- 
tors of SU{2) spin and charge symmetries, respectively. 
In other words, [A^j, 512(21) (i) ^ 512(21) (j)] = q = 

[i?|(^)y,5i"(20)(i)+5i°(20)(j)]. 

Since the su{N) languages provide a complete set of HLs, 
any model Hamiltonian can be written in a similar fashion 
once we identify the appropriate HL and apply the corre- 
sponding SW mapping in the fundamental representation (the 
ordering operators will, of course, have a different meaning 
and algebraic expressions). The key point is that the Hamil- 
tonian operator in the HL becomes quadratic in the symmetry 
generators of the hierarchical group. 

The idea behind any MF approximation is to disentangle 
interaction terms into quadratic ones replacing some of the 
elementary mode operators by their mean value. The crux 
of our HMFT is that the approximation is done in the HL 
where all possible local OPs are treated on an equal footing 
and the number of operators replaced by their mean value is 
minimized since the Hamiltonian is quadratic in the symmetry 
generators. In this way, the information required is minimal. 
In mathematical terms, given Oj- = {Ol) Oy + O}- (Oy ) - 

{Ol.){0,.) + {O}. - {Ol.)){0,. - (Oy)), for an arbitrary bond- 
operator Oy, the approximation amounts to neglect the latter 

fluctuations, i.e., Ol-O,^ « (O^j)Oy + 0^ (Oy) - 
[[12]. An important result is that all local OPs are equally 
treated and, moreover, symmetries of the original Hamiltonian 
related to the OPs are not broken explicitly in certain limits. 
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Now we follow Colpa | |15|] and diagonalize para-unitarily the 
Hamiltonian matrix H [p^6p. The application of a homoge- 
neous linear transformation leads to rpll 



H 



5 

keRBZ j=0 



(7) 



where the mode energies tj^k are, at least, two-fold degen- 
erate Jit]]. In Fig. |l[ we display the orders A and B as 
a function of po = i^Ej("jo) at very low temperature 
(/3 = 10) and different ratios of the competing interac- 
tions J/t for a two-dimensional lattice 
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between the OPs of the original problem, Eq. (|l]), and A 
sJ)mf oc A^, and justifies the la- 



1 



and B is given by 
^^..e^Q.(n-r,)^,+ ,. ^^^^ 

beling of the phases AF (antiferromagnet) and BEC (Bose- 
Einstein condensate) in Fig. [l|. A way to qualitatively under- 



FIG. 1 : Order parameters as a function of the density po for different 
values of J/t and inverse temperature (3 — 10 (in units of t^^). The 
filled circle on the density axis indicates a quantum critical point. 

The resulting MF Hamiltonian together with the SW-boson 
constraint (with Lagrange multiplier X) H = H]\jp + 
^ + "jo + ?^jT) reads [|l 

(ij) '^=hl 

^ njo + A ^(nji + njo + nji) 
j j 

= E I^^^kT^-k+Qi+As (6|,-rfoLk+Qo + 

keRBZ 

^ki^-k+Qo) + H-c] + (A - /x)nko +A^ rik,^ (6) 

where the sum of momenta k is performed over the re- 
duced Brillouin zone (RBZ) with AF ordering wave vector 
Q and rtka = ^ka^ka' ^i*-^ ^ka representing Fourier trans- 
formed modes. = —2jA'yi^, Ab = — 4ti?7k, with 
7k = - e^^'^ (6 are nearest-neighbor vectors). Note that 
when B = in Hmf, the SU{2) spin and U{1), S°°, sym- 
metries are conserved; the opposite case A = preserves 
510(01) +^20(02) 511 ^^22 _500 symmetries. In Eq. (|) 

we have only considered homogeneous solutions [0. 

The corresponding self-consistent MF equations to solve 
are 

keRBZ 
g ^ 

< ^ = ^ E Tk(6L^-k+Q0>A/i=^ ' 
* keRBZ 

" keRBZ a 

We are thus left with a non-interacting system of SW bosons. 



tOis7i(2) 




FIG. 2: Schematics of the order parameter change (as a function of 
Po) with the three su{2) axis describing 3 different su{2) subalge- 
bras of sit(3). Note that the plane j — J, represents the 5 charge 
symmetry generators while the || axis is associated to the remaining 
three generators of magnetism. 

stand this quantum phase diagram is to look into the OF space 
as displayed in Fig. ^. As po (or chemical potential) varies 
from to 1, the OP (depicted as an arrow) moves in OP space. 
When Po = the order is purely AF and the arrow lies on 
the t j axis. For po 7^ 0, the arrow has projections onto the 3 
su{2) axis, i.e., the AF state coexists with a BEC state. There 
is a particular critical value of po ~ poc < 1 for which the AF 
ordering vanishes and the OP is purely BEC with the arrow 
lying in the t — | plane. This BEC ordering persists until 
Po = 1, where it vanishes. 

As can be inferred from our presentation, there are two 
complementary aspects to studying competition and coexis- 
tence between phase orderings in strongly coupled quantum 
systems. One is the direct discovery of the hidden unity 
and subsequent determination of the possible phases and their 
transitions, given a Hamiltonian operator modeling the com- 
plex material of interest. This is the aspect we have described 
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in the present paper The second aspect, to be discussed in 
a separate pubUcation, involves the design or engineering of 
new states of matter using the inverse path of logic. Essen- 
tially, the idea consists of tailoring effective Hamiltonians 
based upon a general symmetry analysis of the possible or- 
derings one would like to realize at zero temperature. Tun- 
ing the parameters of these symmetry-based effective Hamil- 
tonians allows one to move in parameter space along the pre- 
viously established orderings. Indeed, this strategy finds its 
experimental realization in recent work done on atomic BEC 
systems in optical lattices [^, and our approach provides a 
unique theoretical guidance to achieve that goal. 

There are some open issues. One regards the application of 
the HMFT approach to study fermionic problems, for exam- 
ple, a Hamiltonian like Eq. but where the operators for 
different modes on a lattice i anticommute. The method could 
certainly be used, however, fermions do introduce a non-local 
gauge potential [ [To| ] leading to an effective dynamical frus- 
tration which is difficult to handle in a controlled manner 
Another issue concerns the application of the HMFT method 
when longer-range interactions are involved. There is already 
evidence from work on the J1-J2 SU{2) Heisenberg model 
[ [T^ ] that our HMFT will work in those cases as far as homo- 
geneous phases are concerned. Actually, the SW MF theory 
introduced by Arovas and Auerbach [ poj ] in the fundamen- 
tal representation is a particular case of our general HMFT. 
Finally, our methodology exhausts all broken symmetry in- 



stances but it is still quite possible to have purely topologi- 
cal quantum orders and their corresponding phase transitions 
which cannot be described by broken symmetries and associ- 
ated OPs [|T| and, thus, are not included in our framework. 

Summarizing, we developed a theoretical framework and 
a calculational scheme to study coexistence and competition 
of thermodynamic phases in strongly correlated matter In 
our method (given a Hamiltonian modeling the physical sys- 
tem) the order parameters are not guessed but rigorously de- 
termined from group theoretical considerations as symmetry 
generators of a hierarchical language. In this way, the Hamil- 
tonian operator (which does not necessarily have the full sym- 
metry of the hierarchical group) is expressed in terms of sym- 
metry generators. Then, in a non-phenomenological approach 
dubbed hierarchical mean-field theory, we approximated the 
dynamics (and thermodynamics) treating all possible local or- 
der parameters on an equal footing, i.e., without preferred 
symmetry axis. One could say that this procedure follows 
the guiding principles of maximum symmetry and minimum 
information. This allowed us to obtain in a simple manner 
the phase diagram of a model problem exhibiting coexistence 
and competition between antiferromagnetism and superfluid- 
ity. Combined with an analysis of the fluctuations (to analyze 
the stability of the mean-field) one now has a simple machin- 
ery to design phase diagrams. 

This work was sponsored by the US DOE. We thank J.E. 
Gubernatis for a careful reading of the manuscript. 
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